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Abstract 

Extending the idea formulated in Makino et a/[Phys.Rev.E 67,066205], that is 
based on the Berry-Robnik approach [M.V. Berry and M. Robnik, J. Phys. A 17, 
2413], we investigate the statistical properties of a two-point spectral correlation for 
a classically integrable quantum system. The eigenenergy sequence of this system is 
regarded as a superposition of infinitely many independent components in the semi- 
classical limit. We derive the level number variance (LNV) in the limit of infinitely 
many components and discuss its deviations from Poisson statistics. The slope of the 
limiting LNV is found to be larger than that of Poisson statistics when the individual 
components have a certain accumulation. This property agrees with the result from the 
semiclassical periodic-orbit theory that is applied to a system with degenerate torus 
actions[D. Biswas, M.Azam,and S.V.Lawande, Phys. Rev. A 43, 5694]. 
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1. Introduction 



Study of energy level statistics has played an important role in elucidating the universal 
properties of quantum systems, which in the semiclassical limit, reflect regular and chaotic 
features of classical dynamics. Berry and Tabor conjectured that the eigenenergy levels of a 
quantum system, whose classical dynamical system is integrable, have the same fluctuation 
property as uncorrelated random numbers from a Poisson process, and thus are character- 
ized by Poisson statistics^ This conjecture is in contrast with the conjecture of Bohigas, 
Giannoni, and Schmit, which assert that GOE or GUE statistics based on random matrix 
theory (RMT) are applicable to the fluctuation property of energy levels of a quantum sys- 
tem whose classical dynamical system is fully chaoticP These contrasting conjectures have 
been examined using various statistical observables, e.g., the nearest-neighbor level-spacing 
distribution (NNLSD) and the level number variance (LNV).'^^ 

The NNLSD P{S) is the observable most commonly used to study short-range fluctua- 
tions in a spectrum. For eigenenergy levels on the unfolded scale,^ this observable is defined 
as the probability density of finding a distance 5* between adjacent levels. In Poisson statis- 
tics, it is characterized by the exponential distribution P{S) = e~^] on the other hand, in 
GOE/GUE statistics, it approximates the Wigner distribution very well.'^ 

The LNV S'^{L) is the observable commonly used to study correlations between pairs of 
levels, which characterizes long-range fluctuations in a spectrum. It is defined as the average 
variance of the number of levels in an energy interval containing an average of L levels. On 
the unfolded scale, this interval is equivalent to an interval of length L, and the LNV is thus 
defined as 

IJ\L) = me + L)-Af{e)-Lf), (M) 

where A/'(e) is the number of eigenenergy levels below e, and the brackets (■ ■ ■ ) represent 
the average over the value of e. In Poisson statistics, the LNV is equal to the number it- 
self [i.e., IJp^.^^^^{L) = L]; on the other hand, in RMT, in which /? = 1 and 2 correspond to 
the GOE and GUE level statistics, respectively, the LNV increases logarithmically with L 
[i.e., S^iL) ~ (2//?7r^) ln(27rL)] According to Berry's semiclassical theory,'^ the LNV of 
a quantum system with / degrees-of-freedom should display these universal statistics in the 
region L « L^^^ ~ h^^-f~^^ as a consequence of the uniform distributed property of periodic 
orbits in the phase space.® In the semiclassical limit, in which the Planck constant tends to 
zero {h ^ 0), one expects to observe the universalities for all L > in systems with / > 2. 

There are many works examining the Berry- Tabor conjecture in classically integrable 
quantum sjstems^^^^^^^^^^^^^^ Although the mechanism that sup- 
ports this conjecture remains to be clarified; the statistical property of eigenenergy levels 
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to be characterized by Poisson statistics is now widely accepted as a universal property of 
generic integrable quantum systems in the semiclassical limit. 

One possible mechanism producing Poisson statistics has been proposed by Makino et 
al^ based on the Berry- Robnik approach.!^ We briefly review the outline below: For an in- 
tegrable system, individual orbits are confined in each inherent torus whose surface is defined 
by holding its action variable constant, and the whole region of the phase space is densely 
covered with invariant tori. In other words, the phase space consists of infinitely many re- 
gions, which have infinitesimal volumes in Liouville measure. Because of the suppression of 
quantum tunneling in the semiclassical limit, the Wigner functions of quantal eigenstates are 
expected to be localized in the phase space region explored by a typical trajectory,'23''23 and 
to form independent components. For a classically integrable quantum system, the Wigner 
function localizes on the infinitesimal region in — > and tends to a 5 function on a torus 
that is designated by quantum numbers.'^ Then, the eigenenergy levels can be represented 
as a statistically independent superposition of infinitely many components, each of which 
gives an infinitesimal contribution to the level statistics. Therefore, if the individual spectral 
components are sparse enough, one would expect Poisson statistics to apply here, as a result 
of the law of small numbers.'^ 

The statistical independence of spectral components is assumed to be justified by the 
principle of uniform semiclassical condensation of eigenstates in the phase space and by the 
lack of their mutual overlap, and thus, can only be expected in the semiclassical limit.'^''^ 
This mechanism was introduced as a basis for the Berry-Robnik approach to investigate 
the relation between the phase space geometry and the formation of energy level sequence 
in a generic mixed quantum system, whose classical dynamics is partly regular and partly 
chaotic.'^ The validity of the Berry-Robnik approach is confirmed by numerical computa- 
tions for the mixed quantum systems in the extremely deep semiclassical region,'^''^ which 
is called the Berry-Robnik regime^ 

Based on this view, Makino and Tasaki investigated the short-range spectral statistics 
of classically integrable systems.'^ They derived the cumulative function of NNLSD, i.e., 
M[S) = P{x)dx, is derived in the limit of infinitely many components, which is charac- 
terized by a single monotonically increasing function /i(0; S) of the nearest level spacing S 
as 

M^(5) = l-[l-/i(0;5)]exp(^-^ [1 - /i(0; a;)] rfx^ . (1-2) 

The function /2(0; 5") classifies Mp{S) into three cases: Case 1, Poisson distribution Mp{S) = 
1 — for ^S* > if /i(0; +oo) = 0; Case 2, asymptotic Poisson distribution, which con- 
verges to the Poisson distribution for S +oo, but possibly not for small spacings S if 
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< /i(0; +00) < 1; and Case 3, sub-Poisson distribution, which converges to 1 for S* +00 
more slowly than does the Poisson distribution, if /i(0; +00) = 1. Therefore, the Berry- 
Robnik approach when applied to classically integrable quantum systems, admits deviations 
from Poisson statistics. 

Cases 2 and 3 are possible when the eigenenergy levels of individual components show 
a singular NNLSD result from strong accumulatioiJ^ (see also Section [3] of the present 
paper), which is expected when the classical dynamical systems have a spatial symme- 



lj.^mmmmmmmm or a tlme-reversal symmetryP'IS'EB One possible example 



is a rectangular billiard with a rational ratio of squared sidesP''^''^''^''^ These results 
suggest the existence of a new statistical law different from Poisson statistics in the strongly- 
degenerate quantum systems, and raises the question as to whether similar behaviors appear 
in the two-point spectral correlation. 

Deviation from Poisson statistics, resulting from the symmetry, is also exhibited within 
the framework of the periodic-orbit theory.'^S Based on the semiclassical theory of Berry 
and Tabor,'^ Biswas, Azam, and Lawande investigated the LNV for classically integrable 
systems with degeneracy in orbit actions.l^ They showed that the slope g^-^{L) of the LNV 
Z'^(L) = g^^{L)L is described by the average degeneracy of actions of the periodic orbits, 
which for the degenerated quantum systems, is greater than 1, and it is the slope of the 
variance of Poisson statistics. This result leads to the possibility that the degeneracy of ac- 
tions, induced by the symmetry, could be an essential factor for the singularity of individual 
spectral components that yields Cases 2 and 3. This possibility is confirmed by examining 
the properties of LNV for Cases 2 and 3. 

In this study, extending the theory of Makino and Tasaki,'!^ we investigate the LNV 
Z'^(L) of quantum systems whose energy levels consist of infinitely many independent com- 
ponents, and show that the non-Poisson limits (Cases 2 and 3) are possibly observed also in 
the long-range spectral fluctuations. 

Based on the Berry-Robnik approach, the overall LNV is derived as follows: We con- 
sider a system whose classical phase space is decomposed into disjoint regions that give 
the distinct spectral components. The Liouville measures of these regions are denoted by 
Pn{n = 1, 2, 3, ■ ■ ■ , A^), which satisfy X]n=i Pn = 1- Let E{k; L), k = 0,1,2, ■■ ■ be the distri- 
bution function, which denotes the probability to find k levels in an interval (0, L).'2^'E3S 
The LNV S'^{L) is expressed by E{k; L) as 




(1-3) 



fc=0 
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Let P{k; S), k = 0,1,2, he the level-spacing distribution, which denotes the probability 
density to find k levels in an interval of length S beginning at an arbitrary level and ending 
at the level ej+fc+i- P{k; S) is related to E{k; L) as 

P{k; L) = ^ Y.^k - J + l)E{r, L), (14) 

3=0 

and 

/ + 00 f + OO 

dx J [Pik;S)-2P{k-l;S) + Pik-2;S)]dS, (1-5) 

where P{j < 0; S) = oP'E^'EHJ' Eqs. fll-4p and (II -51) are known as the formulae in the the- 
ory of point process, and are derived as corollaries of the Palm-Khinchin formula.'^ The 
NNLSD, which was expressed by P{S) at the beginning of this section, is the special case 
with k = 0. 

When the entire sequence of energy levels is a product of statistically independent su- 
perposition of N subsequences, E{k;L) is decomposed into those of subsequences, en{k; L) 
a^S 

N 

Ej^{k;L)= ^ JJe„(A;„;p„L), (1-6) 

n=l ''I— « 

where e„ satisfies the normalization conditions 

+ 00 

^e„(A;;p„L) = l (1-7) 

fc=0 

and 

^ ken{k; PnL) = pnL. (1-8) 

fc=0 

In terms of normalized level-spacing distribution pn{k; S) of the subsequence, en{k; L) is 
specified as 



en{k;L)=pnJ dx J [pn{k; S) - 2pn{k - 1; S) + Pn{k - 2; S)]dS, (1-9) 

where Pn{j < 0; S) = 0, and Pn{k; S) satisfies the normalizations 

"+00 

Pn{k;S)dS = l (MO) 



and 

k + 1 



Spnik;S)dS = (Ml) 

Pn 
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Note that in general, the individual components are not always unfolded automatically even 
when the overall spectrum is unfolded. However, in a sufficiently small interval [e, e + Z\e], 
each spectral component obeys the same scaling law (see Appendix of Makino et ai^^') and 
thus is unfolded automatically by an overall unfolding procedure. In the Berry-Robnik ap- 
proach. Eg. (11-61) relates the level statistics in the semiclassical limit with the phase space 
geometry. 

In most general cases, the individual components might have degeneracy of levels that 
lead to singular level spacing distributions. In such a case, it is convenient to use its cumu- 
lative distribution function /i„(A;; S), 

Hnik;S)= / pn{k]x)dx. (1-12) 



This function satisfies 

^n{k — 1; 5*) > Unik] S) for all k >1 and for all 5* > 0, (1'13) 
and the normalization condition 

+ 00 

W{k - 1; 5) - /x„(fc; S)] = 0, (1-14) 

fc=0 

with /x,„(j = -1; S) = 0, since /„(fc; S), 
fn{k;S) = 



l^n{k - I; S) - jjnik; S), k>l 
l-fin{0;S), k = 



denotes the probability to find k levels in an interval of length S beginning at an arbitrary 
level, and obviously satisfies, from the definition of this function, the conditions fnik; S") > 
and T.k=o fn{k; S) = 1. 

In addition to Eg. (11-61) . we introduce the following two assumptions that were introduced 
in Ref.P 

Assumption (i). The statistical weights of independent regions vanish uniformly in the 
limit of infinitely many regions: 

maxp„— »-0 as — ^ +oo. (1'15) 

n 

Assumption (ii). For A; = 0, 1, 2, ■ ■ ■ , the weighted mean of the cumulative distribution of 
energy spacing, namely, 

N 

fi{k; x) = Y PnfJ'nik; x) (1-16) 

n=l 
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converges as ^ +00 to ^{k; x), 



lim fi{k;x)=fi{k;x), 



(1-17) 



where the convergence is uniform on each closed interval: < x < S*. In the Berry-Robnik 
approach, the statistical weights of individual components are the phase volumes (Liouville 
measures) of the corresponding invariant regions. 

Under the Assumptions (i) and (ii), Eqs. fll-3p and fll-6p . we obtain the overall LNV in 
the limit — ^ +00: 



When S) = for all k, the LNV of the whole energy sequence reduces to = 
L. This condition is expected when the individual components are sufficiently sparse. In 
general, we expect fi{k; S) > 0, which corresponds to a certain accumulation of the levels of 
individual components. In this case, the LNV of the whole energy sequence deviates from 
^poisson(-^) such a way that the slope of the variance is greater than 1. 

The present paper is organized as follows. The limiting LNV fll-lSp is derived from 
Eqs. (ll-3p and (II -Op . and Assumptions (i) and (ii) in Section [2l In Section [3l the property of 
the limiting LNV is analyzed for the Cases 1-3, wherein the deviations from Poisson statistics 
are clearly observed in Cases 2 and 3. We present a numerical analysis for the rectangular 
billiard that shows deviations from Poisson statistics in Section HI In the concluding section, 
we discuss some relationships between our results and other related works. 



In this section, by using Eqs. fll-61) and (11 •31) . and Assumptions (i) and (ii) introduced in 
SectioEdl we derive the limiting LNV, 



in the limit of infinitely many components A^ — ^ +00. 

First, we rewrite Eq. (ll-3p in terms of the function e„, and decompose it into the LNV 
of individual components: 




(M8) 



2. limiting level number variance 




(2-1) 




(2-2) 



N +00 



N 




n=l k=0 



n=l 
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where we have used the properties Xlfc^ E^ik] L) = 1, Xlfc^ kEj^ik] L) = L, and 

+00 N +00 N 

J2 k'E^ik; L) = J2Y1 ^nenikn, PnL) + Y.{PnL)\ (2-4) 

fc=0 n=l k„=0 n=l 

which follow from Eqs. fll-71) and (11 -SI) , and the relationships k = "^^=1 '^n=i Pn = ^ 

(see also Appendix A). The functions e„(/c; pnL) in the above equation are rewritten in terms 
of the cumulative level-spacing distribution functions Pn{k', S) of individual components. 

'Pn [Pn{k; S) - 2fin{k -1;S)+ fin{k -2;S)], k>2 

e^ik; PnL) = \ -p„ dS [/i„(l; S) - 2/i„(0; 5) + 1] , k=l (2-5) 

Pn dS [I - piri{^] S)] , k = 

Pn dS [pn{k; S) - 2pn{k -1;S)+ pn{k - 2; S)] , k>2 
Pn dS [/i„(l; S) - 2/i„(0; 5) + 1] , k = l (2-6) 

- Pn j'l^ dS [I - PnO-] S)] , k = 

where Eq. (12-61) follows from Eq. (ll-lip . integration by parts and the following limits which 
result from the existence of the average: 

lim S [pn{k; S) - 2p{k - I; S) + pn{k - 2; 5)] = 0, k>2 

lim S [pn{l] S) - 2p{0] S) + 1] = 0, k = l (2-7) 

lim S [1 - pn{0; S)] = 0, k = 

Then, a^{pnL) is described in terms of pn{k] S) as 

+00 

(^lipnL) = k^Cnik; PnL) - pIL^ 

k=0 

= Pn dSy2i2k + l)[pn{k-l;S)-pnik;S)]+Pn [l - ^„(0; ^)] rf^ - p^L^ 
Jo ^^-^ Jo 

= PnL + 2 dSY,PnPn{k-S)-plL\ (2-8) 

where we have used Eq. (ll-14l) and the relation 

+00 

[{k - l)Pn{k - 1; ^) - kpn{k- S)] = (2-9) 

A:=0 

with /i„(— 1; S) = 0. Therefore, in the limit — > +00, we have the convergence 

/.L +00 



r|(L)= hm Y^l(PnL)=L + 2 dS^mS), 



(2-10) 
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where the hmit (12- 101) follows from Assumption (ii) and the following property results from 
Assumption (i): 

N N 

EpI < max p„ p„ = max p„ ^ 0. (2-11) 

n=l n=l 

We note the spectral rigidity A^lL), which was introduced by Dyson and Mehta.'^ 
In combination with the LNV, this quantity has played a major role in the study of the 
long-range spectral statistics. According to Pandej^ and Eg. (12- 101) . the spectral rigidity is 
described in terms of Ji{k; S) as 

-L 



= — + — / dxiL^ -2L^x + x^) dSj2j^ik;S). (2-12) 

§3. Properties of limiting level number variance 

We analyze the slope of the limiting LNV 

g[L) = -El{L) = l + - dSY^mS)^ (3-1) 

fc=0 

which has the following convergence: 

+ 00 

lim ^(L) = 1 + 2 V/i(A;;+oo) > l + 2/i(0;+oo). (3-2) 



L— »+oo ' 

fc=0 



Since fin{k', S) is monotonically increasing for S" > 0, monotonically decreasing for k = 
0, 1,2,---, and < Unik; S) < 1, fi{k; S) has the same properties and is bounded by 
/2(0; +oo) as 

< ^) < /i(0; ^) < /i(0; +CX)) < 1, (3-3) 

where /i(0; +oo) classifies the cumulative NNLSD fll-2p into three cases: Case 1, the Poisson 
distribution if /2(0; +oo) = 0; Case 2, the asymptotic Poisson distribution if < /i(0; +oo) < 
1; and Case 3, the sub-Poisson distribution if fl{0; +00) = 1. 

Then, the property of the limiting LNV is evaluated for Cases 1 — 3 as follows: 
Case 1, /i(0; +00) = 0: The limiting LNV ^^{L) agrees with the LNV of Poisson statis- 
tics: Sp^.^^^^{L) = L. Note that this condition is equivalent to "S") = since j2{k; S) 
is monotonically increasing for 5* and decreasing for k. 

Case 2, < /i(0; +00) < 1: possibly deviates from IJp^.^^^^^{L) in such a way that 

the slope of the limiting LNV is 1 at L = 0, increases monotonically with L, and approaches 
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a number 1 + 2/2(0; +00) or more as L ^ +00. 

Case 3, /i(0; +00) = 1: possibly deviates from Ep^.^^^^{L) in such a way that the 

slope of the limiting LNV is 1 at L = 0, increases monotonically with L, and approaches a 
number 3 or more as L — > +00. 

One has Case 1 if the NNLSD of individual components are derived from the scaled 
distribution functions ipn{0; S) as 

HniO;S)=Pn[ iPniO; PnX)dx, (3-4) 

where V5„.(0; PnS) = Pn(0; S)/pn and satisfy 







DO /■ + 00 

(Pn{0; x)dx = 1, / xipniO] x)dx = 1, (3-5) 
Jo 

and are uniformly bounded by a positive constant D : |v5n(0; S)\ < D { 1 < n < N ) . Indeed, 

the following holds: 

Af „5 N N 

S)\ <'^pI / 1^9^(0; pnx) \ dx < PS y^Pn < PS maxpn 

71=1 ^ n=l n=l 

(3-6) 

Such a bounded condition is possible when the individual spectral components are sparse 
enough. 

In general, one may expect Cases 2 or 3 with p(0; S) > 0, which corresponds to strong 
accumulation of energy levels, leading to a singular NNLSD of the individual components. 
Such accumulation is expected to arise from the symmetry of the system. In the next section, 
we will analyze the LNV of the rectangular billiard systems which is known to deviate from 
Poisson statistics. 

§4. Rectangular billiard system 

We present our results on the various statistical measures discussed in the previous section 
for a rectangular billiard system whose spectral statistics has been precisely analyzed in a 
number of works .'^''^''^''^''^''^''^''^ The eigenenergy levels of this system are given by 

Cn/m = n^ + arn^, (4-1) 

where n and m are positive integers, and a is denoted by the lengths of two sides a and b 
as a = a? jh^. The unfolding transformation {en,m} {cn.m} is carried out by using the 
leading Weyl term of the integrated density of states, A/"(e), as 

vr 

^n,Tn (^n,m) l=(^n,m- (^'2) 

AW a 
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Berry and Tabor observed that the NNLSD of this system possibly deviates from Poisson 
statistics when a is rational.^'' In this paper, we study irrational cases in addition to a 
rational case (a = 1) that are described by a finite continued fraction of the golden mean 
(v^+l)/2, 

a = l + — ^ i ^ = ,1,1 + 5], (4-3) 

1+1+ l+l + () ^ ^ ^ ^ 

with an irrational truncation parameter b G [0, 1). 

Fig.l shows the plots of the LNV S'^{L) for a corresponding to the (a) 25th, (b) 8th, 
and (c) 4th approximations of the golden mean, and (d) a = 1. Our analysis is valid in 
the region L << where L^^^^ = A/vrinm«~^^^ for the rectangular billiardP'^^ We 

used energy levels en,m ^ [4000 x 10'^, 4001 x 10'^], which correspond to L^^^ ~ 3.1 x lO''. 
The numerical computation was carried out using a double precision real number operation. 
When the continued fraction is close to the golden mean, S'^{L) is well approximated by 
^poisson(-^) — L [plot (a)]. On the other hand, in cases in which the continued fractions are 
far from the golden mean, 2J'^{L) clearly deviates from Z'p^.^^^^(L) [plots (b) - (d)]. 

Fig. 2 shows the slope g{L) of the LNV for the four values of a corresponding to the plots 
(a)-(d) in Fig.l. When S^{L) is well approximated by Sp^.^^^^^^L) , the slope g{L) is 1 [plot 
(a)]. Since ^k=Qp'i^'i "S") = is equivalent to /i(0; S) = 0, this result corresponds to Case 1 
given in the previous section(see also Eq. fl3-ip ). In the case in which I^'^{L) deviates from 
^poissoiX-^)' slope g{L) is 1 at L = and increases monotonically with L [plots (b) - (d)]. 
In the limit of L — +oo, plot (b) approaches a number less than 3 and this result clearly 
corresponds to Case 2, while plots (c) and (d) approach numbers greater than 3 and these 
results correspond to Case 2 or 3. 

In order to clarify the attribute of plots (c) and (d), we consider 

1 - M{S) 
1 - //(I - M{x))dx 



f^is) = 1 - , (44) 



that is obtained by the cumulative NNLSD M{S) = P{0;x)dx. The function j2{S) is 
equivalent to /i(0; S) in the semiclassical limit e —>■ +00. 

Figs.3(a)-3(d) show — In [1 — M{S)] for the four values of a corresponding to the plots 
(a)-(d) in Figs.l and 2, respectively. The dotted line in each figure corresponds to the 
cumulative Poisson distribution Mp^^^^^^^^S) = 1 — exp(— 5). Plots (a)-(d) in Fig. 4 show 
/i(S') for the four values of a corresponding to Figs.3(a)-3(d), respectively. When S'^{L) 
approximates IJp^.^^^^{L), M{S) fits Mp^i^^^^-^{S) very well[Fig.3(a)]. In this case, p,{S) is 
obviously [plot (a) in Fig.4]. In case that S'^{L) deviates from Up^^^^^^{L), M{S) for small 
value of S clearly deviates from Mpoiggo„(S') [Figs.3(b)-3(d)]. However, for large value of S, it 
approaches a line whose slope is l(see the dashed line in Figs. 3(b)-3(d)). In these cases, 
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fi{S) in S* ^ +00 approaches a number /i(+oo) such that < /i(+oo) < 1 [plots (b) - (d) 
in Fig.4]. Therefore, plots (b)-(d) correspond to Case 2. 

For the rectangular billiard in the finite energy region, we have not yet succeeded in 
observing Case 3. This case is expected to arise in a square billiard (a = 1) in the high 
energy limit where a stronger accumulation of levels is generated. Based on the number- 
theoretical result of Landau (1908),S3 Connors and Keating have proved that the eigenenergy 
levels en,m of square billiard show a logarithmic increase in the mean degeneracy of levels as 
e — > +CXD, which is described as 

1 -M(+0) ~ --^ ^ 0, (4-5) 
TT In e 

where C2 converges to give C2 — 0.764.'^ The above limit corresponds to a delta function 
of NNLSD, P{S) = 5{S), and is consistent with /i(0; S = +0) = lime_+oo M(+0) = 1, which 
indicates an extremely slow approach to Case 3 in e ^ +00. 

Fig. 5 shows the cumulative NNLSD M{S) for a = 1. This function is not smooth at the 
level spacings separated by a step tt/A^ which correspond to accumulation of levels. Note 
that M(+0) > due to the degeneracy at 5* = 0. Since /i(+0) = M(+0) > 0, this degeneracy 
at S" = is identified also in Fig.4. As the eigenenergy levels become higher, M(+0) increases 
monotonically and approaches 1. In the limit e —>■ +00 where M{+0) = /2(0, 5* = +0) = 1, 
all steps except the step at S* = are suppressed since M{S) is monotonically increasing for 
^> 0. 

Fig.6 shows 1 - M(+0) vs AC2/iT^/\ne for various energy ranges. Although we are 
not yet far enough in the high energy region where 1 — M(+0) = 1 — /i(+0) << 1, the 
agreement between them is very good, and is better as e — +00. Therefore, the extremely 
slow convergence to Case 3 is well reproduced by a numerical computation. The almost same 
results for a = 22/21 have already been reported by Robnik and Veble.'® 

§5. Conclusion and Discussion 

Based on the approach of Berry and Robnik, we have investigated the energy level statis- 
tics of classically integrable quantum system and discussed its deviations from Poisson statis- 
tics. In the Berry-Robnik approach, individual eigenstates localizing on the different phase 
space regions provide mutually independent contributions to the statistics of energy levels 
in the semi classical limit. Since the phase space of integrable system is densely covered with 
invariant tori, the eigenfunctions of the classically integrable quantum system are localized 
on the regions that have infinitesimal volumes in Liouville measure. Therefore, we have 
considered the situation in which the eigenenergy sequence is a superposition of infinitely 
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many independent components, and each of which gives an infinitesimal contribution to the 
level statistics. Moreover, by developing the approach of Makino et al^ into the statistics 
of higher order (long-range) spectral fiuctuations, the LNV of systems consisting of infinitely 
many spectral components are obtained. The LNV is characterized by the monotonically 
increasing functions fi{k; S),k = 0, 1, 2, ■ ■ ■ , of the level-spacing S, where the lowest order 
term /2(0; 5*) is associated with the NNLSD. The property of the LNV is classified as follows: 
Case 1, /i(0; +oo) = where the cumulative NNLSD is the Poisson distribution, the LNV is 
the Poissonian Z'p^j^^^^(L) = L; Case 2, < /2(0; +oo) < 1 where the cumulative NNLSD is 
the asymptotic Poisson distribution, the LNV deviates from Z'p^.^^^^^(L) in such a way that 
the slope is greater than 1 and approaches a number > 1 + 2/2(0; +oo) as L — +oo; Case 
3, /2(0; +oo) = 1 where the cumulative NNLSD is the sub-Poisson distribution, the LNV 
deviates from Ep^.^^^^^{L) in such a way that the slope is greater than 1 and approaches a 
number > 3 as L — +oo. Therefore, we have shown that deviations from Poisson statistics 
(Cases 2 and 3) are possibly observed, not only in the property of NNLSD characterizing 
short-range spectral fiuctuation as shown in the work of Makino et al^ but also in the 
property of the LNV characterizing the fiuctuations of all ranges. 

Note that Cases 2 and 3 may arise when there is a strong accumulation of levels, which is 
characterized by the singular NNLSD of individual components. Such accumulation would 
be expected when there is a spatial symmetry or time-reversal symmetry. One example is 
a rectangular billiard shown in Section H] where the result shows Case 2 in addition to Case 
1, and an extremely slow approach to Case 3 . Similar results due to a spatial symmetry 
are also shown in the equilateral-triangular billiard,'^''® torus billiard,'^''® and integrable 
Morse oscillator. Another example is a certain type of systems with the time-reversal 
symmetry studied by Shnirelman,'^ Chirikov and Shepelyansky,'33J and Frahm and Shep- 
elyansky.l^ In this system, the strong accumulation of energy levels, resulting from the time 
reversibility, is refiected in the NNLSD as a sharp Shnirelman peak at small level-spacings. 

Rigorous results confirming the Berry-Tabor conjecture (Case 1) are also reported for a 
certain classically integrable system. In the work of Marklo^^ and Eskin et al^ eigenvalue 
problems are reformulated as lattice point problems, and it is exactly proved under explicit 
diophantine conditions that a two-point spectral correlation exhibits Poisson statistics. 

The periodic-orbit theory, applied for an integrable system with degeneracy in orbit ac- 
tions, gives a similar result. In the work of Biswas et al,^"^- it was shown that the slope g^^{L) 
of the LNV S'^{L) = g^^{L)L represents the average degeneracy of actions of periodic orbits, 
which is greater than 1 for the rectangular billiard with a rational ratio a of squared sides. 
Since this property is qualitatively consistent with the property of the limiting LNV fl2-10p 
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in Cases 2 and 3, one might expect to have the relation: 

2 j-L 

g^^{L) = l + - / Y.Ji{k-S)dS, (5-1) 

fc=0 

which enables us to discuss the properties of the overall NNLSD from the periodic-orbit 
theory. From Eqs. (l5-ll) and fl3-3p . we expect that the Case 1 corresponds to the non- 
degeneracy of actions: g^^{L) = 1, while Cases 2 and 3 correspond to the degeneracy of 
actions: (7av(-^) > 1- Therefore, we have confirmed that Cases 2 and 3 are closely related to 
the degeneracy in orbit actions. At the same time, we have only superficially examined the 
possibility that the degeneracy of actions, induced by the symmetry, could be an essential 
factor for the singularity of individual spectral components that yields Cases 2 and 3. This 
part should be investigated in detail in a future work. 

There is another approach similar to the one presented in this paper, in which the LNV 
is derived in the limit of infinitely many components.'^ The LNV is described by the Dyson 
two-level cluster function Y2 as 



S^iL) =L-2 [ (L - S)Y2iS)dS. (5-2) 
^0 



For spectral superposition, Y2 is described in terms of the cluster function of individual 
components ?/2,n as in Pandey's workpS 



TV 



Y2iS)=J2ply2ApnS). (5-3) 

n=l 

Then, for the overall LNV, one has the convergence 

lim S\L) =L + 2 [ c{S)dS (54) 

^ + 00 Jq 



N 

with 



N 



c{S) = - lim y^Pn y2,n{x)dx, (5-5) 

N~^ + QO ^ ' In. 

n=0 

where c{S) = corresponds to Poisson statistics and c{S) 7^ indicates deviations from 
Poisson statistics. Since ?/2,n is associated with the level-spacing distribution of individual 
components as y2,n{x) = 1 — Ylit=QPn{k]x)^ c{S) is rewritten as 



c{S)= lim y^^pn y^Pnik;x)dx. (5-6) 

n=0 k=0 
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When Pn{k; 5), A; = 0, 1, 2, ■ ■ ■ , are assumed to satisfy the commutation relation J^"^ ^k=oPn{^'i S)dS = 
J2k=o lo''^ Pnik; S)dS, c{S) is described as c{S) = J2k=of''i^'j ^) hmit (15-41) is consistent 
with the hmit fl240|) . 
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Appendix A 



Derivation of Eqs. h2-!^) and ^2-4^ 



We briefly show the detailed process to derive Eqs. fl2-3p and fl2-4l) from Eqs. fll-41) . fll-7p . 



and fll-Sp . First, we rewrite Eg. (11-31) as 




(A-1) 



/c=0 



fc=0 



where we have used the normalization conditions 'Ylik=o ^Nik] L) = 1 and X]fc^ kE^i^k; L) = 



L, which follow from Eqs. (ll-7p and (11-81) . and the relations k = J2n=i X]n=i Pn = 1 



as 




(A-2) 




N 




(A-3) 



n=l 



By using Eqs.(jMl) and (IX^ . PEN{k; L) in the right hand side of Eq. dTTT]) is calcu- 



lated in a similar manner: 




= XI 5Z ( n (ii{h]piL)\^klen{kn]pnL) 
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A'' +00 +00 

+ 

n^i kn=0 ki=0 

N / +00 \ \N-1 

XI I 5Z klenikn, PnL) J I ^ d^niki] PiL) J (A-4) 
n=l \fcn=0 / \fei=0 / 

N +00 \ ^ / +00 



+ 

\n=l kn=0 / n=l \fc„=0 / 

N +00 N 

= J211 ^nenikn, PnL) + - Y,{PnLf . (A-5) 

n=l fe„=0 n=l 

Since the off-diagonal terms i vanisfi in tfie second equality, we have 

+00 N / +00 

k^E^ik; L)-L^^J2[J2 ^n^^nikn, PnL) - (PnLf 
k=0 n=l \fc„=0 

N +00 

= (fcn - pnLfenikn, PnL). (A-6) 

n=l fen=0 

Therefore, the overall LNV of the system with N independent components is described as 

N 

SUL)^J2^l{PnL), (A-7) 

n=l 

where cr^ is the LNV of the spectral component: 



(^niPnL) = ^{k - PnLfenik; p^L). {k-i 



k=0 
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L 

FIG.l Level number variance S'^{L) of the rectangular billiard systems for (a) 25th, (b) 8th, 
(c) 4th approximations of a = (a/S + and (d) a—1. The solid line denotes the LNV 
of Poisson statistics, Z'^^.^^^^(L) = L. We used energy levels e„,^ e [4000 x 10''', 4001 x 10''']. 
The total numbers of levels are (a) 10000457, (b) 10000428, (c) 10000266, and (d) 10000162. 
The truncation parameters 5 are (a) tt x 10"^ (b) 7r/3 x 10"'' and (c) tt x 10~^. 




FIG.2 Slope of the level number variance, g{L) = E'^{L)/L, for (a) 25th, (b) 8th, (c) 4th 
approximations of a = (v^+ and (d) a — 1. The sohd hne, g{L) — 1, corresponds to 
Poisson statistics. 
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FIG.3 Function -ln[l - M{S)] for (a) 25th, (b) 8th, (c) 4th approximations of a = + 
l)/2, and (d) a = 1. The dotted hne corresponds to the cumulative Poisson distribution 
M{S) = 1 - e-^. 
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FIG.4 Parameter function jjL{S) for (a) 25th, (b) 8th, (c) 4th approximations of a = (a/S + 
l)/2, and (d) a — 1. The sohd hne, jjL{S) — 0, corresponds to the Poisson statistics. 
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FIG.5 Cumulative NNLSD M{S) for a = 1. We used 999844 levels e„,^ e [40 x 10^ 41 x 10^] 
for plot (A), 10001872 levels e„,^ e [4000000 x 10^4000001 x 10^] for plot (B). We observe 
M(+0) ~ 0.767 for plot (A) and M(+0) ~ 0.826 for plot (B). The dotted curve is the 
cumulative Poisson distribution M{S) — 1 — 
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FIG. 6 Numerical test of Eq.( 14-5p for a square billiard (a = 1). The solid line exhibits the 
theoretical prediction, which is valid in the semiclassical limit e — > +00. In each plot, we 
used 10^ eigenenergy levels obtained from the unfolded energy range. The error bar exhibits 
the energy range acquiring numerical data. 
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